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Abstract 

We present quantum Monte Carlo results for the Hubbard model on a ladder 
using the re-structuring method which employ eigenstates of two-site system on a 
rung to construct a complete set. 

From technical reasons we concentrate on the case in which the hopping along 
the leg is much less than the hopping on the rung. We observe the ground state 
of the half-filled system is made of singlet-like states described in this paper and 
its first excited state contains one triplet rung. When the system is lightly doped 
states composed of an electron and a hole appear on some rungs. We find very few 
hole-pairing on the rung upon doping up to the quarter-filled case. 
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Section 1 Introduction 

Lately strongly correlated systems of electrons on ladders with two or more 
legs have become quite intriguing from both theoretical and experimental points 
of view. Various models are numerically investigated by means of exact diagonaliza- 
tion [1], 0, |3| and renormalization group approaches ||, ^ ||. Active studies on the 
quantum spin ladder systems, which correspond to the half-filled fermionic ladders 
with infinite Coulomb repulsion, are also stimulative [7j . Among their results, espe- 
cially interesting are spin gaps, ground states as well as excited states of the systems 
and formation of the hole-pair bound state. Presence of spin gaps is reported for 
several spin systems and for fermion systems such as the t-J-t'-J' model on a single 
ladder (two coupled chains) studied in Ref. 0, the t-J model on a single ladder 
|| and on multi-ladders with even number of coupled chains 0, and the two-chain 
Hubbard Model [|J . In Refs. [[I], [| emergence of the hole-pairing upon doping, which 
would indicate the high T c superconductivity, is also claimed. 

Quantum Monte Carlo techniques based on the Suzuki- Trotter formula, which 
prove useful for many spin systems |8], [| , are difficult to apply to the fermion systems 
because of the negative sign problem except for those on a single chain [K| . Usually 
studies on the two-dimensional fermion systems have been done after integrating the 
fermionic degrees of freedom jlT, T2J, 13, [TJ]]. Yet it would be desirable to simulate 
such electron systems keeping fermionic degrees of freedom and vividly observe how 
electrons and holes behave on ladders. Recently dynamical properties of the half- 
filled Hubbard ladder are studied by Scalapino et al. using grand canonical 
quantum Monte Carlo calculations with maximum entropy analytic continuation. 
In this paper we formulate the Hubbard Model on a single ladder, with or without 
doping, applying the re-structuring approach 0. Here we employ a complete set 
composed of sixteen eigenstates of electrons sitting two ends of a rung, including the 
double occupancy (See the Table). Although our simulations are limited to the case 
in which the hopping along each chain is much weaker compared to the hopping 
between sites on a rung, we obtain there Monte Carlo results with satisfactory 
statistical accuracy on the system's energy as well as what we call state distribution 
in this paper, which is defined in section 2, sorted according to values of the third 
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component of the total spin. At half-filling we observe that in the ground state the 
singlet-like state, the ninth eigenstate in the Table, is sitting on all rungs of the 
system and the first excited state has a clear gap because of the triplet eigenstate 
which appears on one rung of the system. Upon doping the lowest lying state we have 
observed contains not a hole-pair eigenstate but some of one-electron eigenstates in 
addition to the singlet-like eigenstate. Here we see no spin gap. All our results are 
compatible with the corresponding ones in Ref. [|J. For the half-filled ladder our 
observations are also in agreement with the t-J results in Ref. @. 

In the next section we sketch the model and the method used in this paper. In 
section 3 the Monte Carlo results are presented. The final section is devoted to the 
summary. 

Section 2 the Model 

The model we study is the Hubbard ladder with two legs and N rungs. Its 
Hamiltonian is 

AT-l 

« - -t c E E [(c\ a (i)c la {i + l) + cl(i)c 2(r {i+l))+R.c. 

i=l <r=U 

N 2 N 

n itW™jiW (i) 

»=i o-=t,J. i=l »=1 

where c^ a (i) is the creation operator for an electron of spin a on the z-th rung of the 
j-th leg of the ladder and rij a (i) = Cj a (i)cj a (i). From technical reasons we employ 
an open boundary condition with a hole at each edge of the legs (See Fig. 1). The 
partition function Z is given by 

Z = tr{e^ n -^}, (2) 

with the inverse temperature /3, the chemical potential /i and 

2 N 

^=EEE^« (3) 

j=i t=i CT=T> j 

Let us describe our re-structuring approach briefly. In our representation the 
states of the system are given by 

I a) =| ai, a 2 , a 3 , a N ), \ (3) =| b u b 2 , b 3 , b N ), 
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where a m (b m ), which is one of the sixteen eigenstates of the two-site system shown 
in the Table, denotes the state on the m-th rung along the chain A(B) in Fig. 1. 
We expect this choice is proper at least for strong hopping on rungs. It should be 
noted that a m and b m may differ by sign for the same occupants on the two sites of 
the m-th rung because the left site on the rung along the chain A is the right site 
along the chain B and vice versa. 

Then we divide the above Hamiltonian ([I]) with even N into four parts, 

N/2 

hi [o — t, 

4 



1 

n A1 = -n R -t c J2 £ [ c U 2k - 1)^(2*0 + h.c. 

4 fc=l<r=T4 
1 N,2 ~ 1 

n A2 = -n R -t c E £ [4(2^)^(2^ + 1) + h.c. 

k = l <T=T,J. 

1 N/2 

H B1 = -H R -t c J2 E [ct(2k-l)c 2<7 (2k)+E.c. 

fc=l<r=T4 
1 iV/2-1 

U B 2 = jH R -t c E E [cL(2fc)c 1<7 (2fc + 1) + H.c. 



fe=l cr=t,J. 



(4) 



where 



A" 



2 2V 



7lj-rUm.ii u 



Wa^-trE E [4(i)c 2CT (i) + H.c.]+[/EE 

i=l o-=T4 j=l i=l 

Note that the Hamiltonian 7^ + Ha2 (Hbi + ^#2) describes the Hubbard chain 
A(B) with the hopping parameters t c and t r /2 and the Coulomb repulsion U/2. 
Applying the Suzuki- Trotter formula we obtain 

Z n _ tr |[ e -/3(WAi-M'A^)/™ e -/3(WA2-M'A^)/™ e -/3(W S i-^'Ar)/n e -/3(W B2 -/ i 'Ar)/njn| 



£• 


■■ E £-•• E 






<"ol 


(a 2 «-l|<A)| (ftn-il 










I "2) ("2 


1 A>> 




e -P{HBX~n'M)ln | | e -/3(W B2 -^)/n 


/32></? 2 | 


"3) 



(a 2n -2 I e -«wW/0/» I c^X^-i I e -/»CWW*)/» | ao )( aQ | /3 2n „ 2 ) 

<ftn-2 I e-Wn-M/n | ft B _ 1 )( / fl^ 1 | fi "««»-W» | ft) (ft | ^ ( 5 ) 

where n is the Trotter number and // = ///4. 
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For later use let us introduce a hermitian operator O defined by 

= y £\a)C a (a\, (6) 

where C a denotes a real number. Since 

(a\ 

= C ao (a Q | e -M HM -M/« | Q;i ) ) 

we obtain 

tr{0 e -« w -M0} 

= lim tr{C r e -/ 3 (WAi-M^)/n e -/3(WA2-M'A^)/™ e -/3(WBi-M'A^)/n e -/3(W S2 -MW)/nnn| 
n^oo 

= i™E- E E 

(<*ol (<*2n-l| (A) I <An-i| 

(a„ | e -K H "-»'W» | ai )( ai | e -WA2-M/n | a2 )( a2 | 
(/3 Q I e -fl*Bi-MW)/» | ft) (ft | e -^(w S2 -^'A0/n i i ^ 



(^2n-2 | e-Wn-M/n | ft^Xft^ | e *-^/" I /5 )(/3o I « ). (7) 

With the partition function Z n in (|5]) we simulate the system on the N x 4n 
checkerboard with the number of rungs N up to 32 and the Trotter number n = 16 
and 24. We obtain one new configuration by all possible local updates and global 



updates in both the spacial and the Trotter directions using techniques in Ref. |L5 
We do not fix the number of electrons on a ladder so that the system would come 
to the thermal equilibrium swiftly. 

Since one is not free from the negative sign problem with the formulations stated 
above, our study is limited to the parameter regions where the cancellation is not 
serious. To check it we measure the r ratio, 



where Z + (ZJ) is the number of configurations with positive (negative) weight. We 
refer to the part that survived the cancellation as the whole configurations. A 
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physical quantity (A) is calculated by 

(A) = (9) 

with contributions from positively (negatively) signed configurations A + (A_). 

The physical quantities we measure are the hole doping 5 and the system's energy 
(E). We obtain 5 by 

A^ e being the number of electrons on the ladder. As for the system's energy, we 
measure the quantity G = E — fiN e using the standard technique in the quantum 
Monte Carlo method and obtain (E) by 

(E) = (G + fiN e ). (11) 

In order to shed some light on the system's wave function we also measure the 
frequency distribution for every eigenstate in the Table, which is the ratio of the 
number of rungs on a ladder occupied by that eigenstate to the total number of 
rungs A" and will be called state distribution. It should be noted this quantity is 
physical too, because the operator to represent the number of each eigenstate on a 
ladder is expressed by (|j). For the hole-pair eigenstate, for example, C a in (|6|) is 
zero if none of a^s in the state | a) are equal to | 00), is one if | a) contains one and 
only one | 00) state among its a^'s, and so on. 

The third component of the total spin of the system, S z , is another measurable 
quantity in this formalism. We will see S z is very helpful to study system's energy 
spectra. 



Section 3 Results 

Now let us show our Monte Carlo results. Typically first ten thousand configura- 
tions are discarded for the thermalization and next a few ten thousand configurations 
are used for the measurement. Throughout the rest of the paper we take t r =l and, 
unless stated otherwise, U=6. 

We concentrate on the case t c < 1, where the negative sign problem is not so 
serious. In this case the ground state at half-filling is expected to be made of singlet 
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rungs. We simulate the system with t c = 0.1, namely t r /t c = 10. In Fig. 2 we plot 
the state distribution on the N = 16 ladder as a function of (3, with values of /i 
adjusted to ensure 5 small enough to realize the half-filled system. Values of the r 
ratio are between 0.58 and 1 for these (3 and \x. We see the mixed state (Ai), the 
ninth state in the Table, dominates as (3 increases to become almost 100% of the 
whole at f3 — 15. We also observe S z — for 99.6% of the whole configurations 
at this value of (3. These results indicate we can treat the state at (3 = 15 as the 
system's ground state. The ratio of the singlet component in the ground state is 2u\ 
(ui being defined in the Table Caption), which is about 0.92 for U — 6 and goes to 
1 as U goes to infinity. We refer to this state as the singlet-like state. 

In order to study the excitation spectrum at half-filling, we examine the value 
of S z at P — 13 and S ~ 0, where we find that 98.85% of the whole configurations 
has S z — and 0.59 (0.56)% has S z — 1 (—1). The state distribution with non-zero 
S z shows that 6.235 (6.244)% of the S z — 1 (—1) part is occupied by the sixth 
(the eleventh) eigenstate, namely by the triplet(l) (the triplet (—1)) state, the rest 
being occupied by the singlet-like state. These results strongly indicate that about 
1.73% of the whole configurations at f3 — 13 (including the triplet(O) state which 
should be about 0.58/98.85 of the S z — part) is in an excited mode in which 
one of the N(= 16) singlet-like rungs in the ground state turns to the triplet one. 
By measuring system's energies of S z — 0, 1 and —1 configurations separately we 
estimate the energy of this excited state, E 1: from energies of S z = ±1 configurations 
and the ground state energy, E , from energy of S z = configurations and E-y. The 
energy gap E\ — Eq thus obtained is 0.608 ± 0.002, which is quite close to the 
difference of eigenvalues — 2/i — (— Ai) ~ 0.606. 

Similar results are obtained on the ladders with different N. Figure 3 plots the 
energy gap as a function of 1/N. The results for t c = 0.1 up to N = 32 indicate 
that the energy gap of this excitation mode, namely the spin gap, remains finite in 
the thermodynamic limit. Results for t c = 0.2 up to iV = 16 and for t c = 0.3 and 
N = 8 plotted together also support the presence of the gap. Measurement of the 
gap for larger ladders or larger t c was unsuccessful because of the decreasing value 
of r. 
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Let us see what happens to the system upon doping, which is done by reducing 
the value of the chemical potential. More cancellation due to the negative sign 
problem is observed at intermediate values of 5, yet it is not difficult to obtain 
statistically meaningful results except for the region near 5 ~ 0.25. Fig. 4 shows 
state distribution at /3 — 13 for various values of S. Here we observe, from half-filling 
(5 = 0) to nearly quarter-filling (5 = 0.5), the singlet-like state diminishes as 5 grows 
and the ratio of the one-electron states increases instead. As is plotted together, only 
the second and the fourth eigenstate, which are even in the site-exchange, appear as 
the one-electron states. Note that for 5 < 0.5 we see very few hole-pairs on rungs. 
For 5 > 0.5 the distribution changes drastically because the number of holes exceeds 
the number of electrons, but we pursue this region no longer. 

Since we do not fix the number of electrons in the updating processes, configu- 
rations with different numbers of electrons can emerge in the simulation even when 
the chemical potential is adjusted to result in integer values of (N e ). In simulations 
at half-filling this causes no problem because we observe that for large values of (3 
very few configurations have N e different from 2N. For the doped system, however, 
we find the contamination is serious at any available f3. In search for the ground 
state of the system with 2(N — N ) electrons (N = 1,2, •••), we therefore have 
to discard configurations with N e ^ 2(N — N ) in the measurements of the state 
distribution and the energy spectra at the cost of statistical accuracy. The results 
thus obtained at (5 = 9 for the TV = 16 ladder with N = 1 (i.e. for 5 = 0.0625) 
show us that the ratio of the second, the fourth and the ninth eigenstates in the 
Table (namely the even one-electron state with up spin, with down spin and the 
singlet-like state) is about 1/16 : 1/16 : 7/8 for S z = whereas it is 1/8 : : 7/8 
(0 : 1/8 : 7/8) for S z — 1 (—1) and contributions from other states are negligible. 
Consistently with these results, we do not see any energy gap here; the energies 
are (E) = -10.69 ± 0.01, -10.67 ± 0.02 and -10.63 ± 0.02 for S z = 0, 1 and -1, 
respectively. It should be noted that these values of (E) are smaller than —10.48, 
the value with t c — directly calculated from the eigenvalues of the rungs. We 
also simulate the N = 32 ladder with N = 2 at (5 = 9 and obtain statistically 
meaningful results on S z — 0, — 1 configurations which suggest that the system has 
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two (one) even one-electron rungs with up spin, two (three) even one-electron rungs 
with down spin and 28 singlet-like rungs for S z = (—1). The energy on this ladder 
is (E) = -21.30 ± 0.05 (-21.29 ± 0.07) for S z = (-1), which indicates that the 
finite-size effect is small. 

Section 4 Summary 

In this paper we carry out Monte Carlo simulations of the Hubbard model on a 
ladder by the re-structuring method we proposed in Ref . |J . Having constructed the 
complete set in the Suzuki-Trotter formula from the eigenstates for two electrons on 
a rung, our method works satisfactorily when the hopping along the chain is much 
weaker compared to the hopping within rungs. At half-filling we obtain numerical 
results which indicate that the ground state is of singlet-like eigenstate only, the 
first excited state contains one triplet rung besides singlet-like ones and the spin 
gap persists in the thermodynamic limit. When the system is doped we observe 
that on some rungs the eigenstates which contain one electron and one hole appear 
instead of the singlet-like ones seen in the half-filled case. 

We did not observe any hole-pair formation in the parameter region we could 
study. Whether the doped system realizes hole-pairing or not on ladders with larger 
hopping parameters along the chain is an open question so far. Although it seems 
difficult to obtain data for such ladders straightforwardly, detailed study on the 
system's wave function which is possible in this formulation would be useful to 
understand properties of the system. 
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Table Caption 

Eigenstates, eigenvalues for the partial Hamiltonian of two-site system on one 
rung, 

2 2 

h- fin e = -t r i c \aC2a + H.c.) + U Y n rt n jl E E 
<r=t,J. j=i o-=T4i=i 

with t r and [/ in the Hamiltonian ([!]) and the chemical potential /i in (|2|), and their 

names used in this study. 

| Si, S r ) in the Table represents the state on a rung whose left (right) site is 
occupied by Si(S r ), where Si(S r ) is either (hole), | (an up-spin electron), J, (a 
down-spin electron) or | (an electron pair). 

Parameters U\ and u<i and eigenvalues Ai and A2 in the ninth and the tenth 
eigenstates are defined as 



Mi 



«2 



\ 



1 + 



u 



u 2 + mi ' 



\ 



u 



U 2 + 16t 2 ' 



Xl = -^-u + ^u 2 + m 2 ), 



X 2 = -(Afi-U- JU 2 + lQt 2 ). 



Figure Captions 
Figure 1 

A schema to show the ladder we study. The open circle on each edge of the ladder 
denotes a hole. The thin solid (dashed) line is to indicate the chain A (B) described 
in section 2. 
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Figure 2 

The state distribution R for the t r — 1, t c — 0.1, C/ = 6 and iV = 16 Hubbard 
ladder as a function of (3 with | 5 \< 7 x 10~ 3 . The Trotter number n is 16. The 
open circles, right triangles, left triangles and down triangles denote the ratios for 
the first, eighth, tenth and sixteenth eigenstate in the Table, respectively. The 
filled diamonds denote the ninth eigenstate which is referred as singlet-like state in 
the paper. The open squares (asterisks) denote the sum of the contributions from 
the second to the fifth (from the twelfth to the fifteenth) eigenstates, which is the 
contribution of all one-electron (three-electron) states. The open diamonds denote 
the contribution of all triplet states, namely the sixth, the seventh and the eleventh 
eigenstates. Statistical errors are within symbols. 
Figure 3 

The spin gap A = E\ — Eq at half-filling with t r = 1 and U=6 as a function ofl/N, 
N being the number of rungs of the ladder. The system's energy defined by ( |TT| ) is 
measured for the ground state which is made of N singlet-like rungs (Eq) and for 
the excited state where (N — 1) rungs are singlet-like and one rung is occupied by 
the triplet (1) (or triplet (—1)) eigenstate (Ex). The open circles denote results for 
t c = 0.1, (3 = 13 and n = 16, the open square for t c = 0.1, j3 = 13 and n = 24, the 
open triangles for t c = 0.2, (5 — 9 and n = 16, the open diamond for t c = 0.3, (3 = 7 
and n = 16. Statistical errors are within symbols except for the t c = 0.3 datum. 
Figure 4 

The state distribution R for the t c = 0.1, U = 6 and iV = 16 Hubbard ladder 
as a function of 5 with (3 = 13. The Trotter number n is 16. The open circles, 
filled diamonds, open squares and open diamonds denote contributions of the hole- 
pair state, the singlet-like state, all one-electron states and all the triplet states, 
respectively. Contribution from the even one-electron states, namely from the second 
and the fourth eigenstates, is also plotted by the crosses. The open triangles are the 
contribution of the rest. Statistical errors are within symbols. 
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Table 



No. 


eigenstate 


eigenvalue 


Name 


1 


1 oo> 





hole-pair 


2 


7i(IT0>+|0T» 


-(// + t r ) 


one electron (up, even) 


3 


7i(IT0)-|0T)) 


-(n - t r ) 


one electron (up,odd) 


4 


^(I10)+|0|)) 




one electron (down, even) 


5 


^(l|0)-|0|)) 




one electron (down,odd) 


6 


ITT) 


-2/i 


triplet (1) 


7 


73(in)+i4T» 


-2(i 


triplet (0) 


8 


^(II0)-|0T)) 


-{2(i -U) 


an electron-pair and a hole 


9 


«i(IU>- 


-UT»+«2(lto)+|ot» 


-Ai 


mixed state (Ai) 


10 


U2(\U)- 


-UT»-tii(llo>+|ol» 


-A 2 


1 i j / \ \ 

mixed state (A 2 J 


11 


III) 


-2/i 


triplet (-1) 


12 


^(im+iTD) 


-(3^-E/"-f r ) 


three electrons (up,even) 


13 


^(IIT)-ITD) 


-(3// - C/ + t r ) 


three electrons (up, odd) 


14 


^(IID+IID) 


-(3// - C/ - f r ) 


three electrons (down, even) 


15 


^(IID-IID) 


-(3(i -U + t r ) 


three electrons (down,odd) 


16 


III) 


-{A(i - 2U) 


electron-pairs 
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